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We present calculations of the thermal and electric linear response in graphene, including disorder 
in the self-consistent t-matrix approximation. For strong impurity scattering, near the unitary limit, 
the formation of a band of impurity states near the Fermi level leads to that Mott's relation holds 
at low temperature. For higher temperatures, there are strong deviations due to the linear density 
of states. The low-temperature thermopower is proportional to the inverse of the impurity potential 
and the inverse of the impurity density. Information about impurity scattering in graphene can be 
extracted from the thermopower, either measured directly, or extracted via Mott's relation from the 
electron-density dependence of the electric conductivity. 



The recent isolation of the two-dimensional one-atom 
thick honeycomb crystal of carbon,— called graphene, has 
generated tremendous interest. Partly because of the po- 
tential of carbon-based nano-scale electronics, but also 
for fundamental reasons^i^^i 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 Re- 
cent experiments 2 * 3 -^ have shown that the charge conduc- 
tance of graphene reaches a universal value 4e 2 //i at low 
temperatures. It has also been found that the conduc- 
tivity is linearly dependent on the electron density. Fur- 
thermore, an unconventional half-integer quantum Hall 
effect has been discovered. 

All these results are nicely in line with theoret- 
ical results based on an effective low-energy theory 
of graphene within which charge carriers are mass- 
less Dirac feTmions, 16 i 17 i 18 i 19 i 20 i 21 . 22 i 23 i 24 i 25 i 26 i 27 The low- 
temperature conductivity of graphene is most probably 
limited by scattering against impurities, such as vacan- 
cies or other imperfections i 21 i 22 ' 23 i 24 ' 25 For strong scat- 
tering, in the unitary limit, a band of impurity states cen- 
tered at the Dirac point is formed i 22 ' 23 i 28 ' 29 i 30 For small 
temperatures, compared with the impurity band width, 
the conductivity is predicte d 19 ' 20 ! 21 ' 22 ! 23 to reach a min- 
imal universal value 4e 2 /(7r/i). Within the same model 
of strong impurity scattering, the conductivity varies lin- 
early with charge density for chemical potential shifts 
that are large compared with the impurity band width. 23 
In contrast, linear dependence is absent for weak impu- 
rity scattering^ 3 - In a slightly different scenario, Coulomb 
randomly distributed scatterers are predicted to lead to 
similar properties. 25 

In this paper we explore further the consequences of 
impurity scattering near the unitary limit, and present 
results for the response to a thermal gradient in addition 
to the more widely studied electric response. In particu- 
lar, we find that the thermopower can provide valuable 
information about impurities in graphene. 

The related Nernst signal, that appears in a magnetic 
field, was studied theoretically in Ref. lU Further prop- 
erties of the thermopower in a magnetic field was stud- 
ied in Ref. [3l| in the strict unitary limit for which the 
thermopower vanishes at charge neutrality. If vacancies 
are strong scattering centers, the impurity potential is 
large, estimated e.g. by graphene's vacuum level. Near 



unitary scattering can occur, although the strict unitary 
limit can in reality not be reached. This has a big im- 
pact on the thermopower, as the impurity band is shifted 
from the Dirac point which leads to a large electron hole 
asymmetr y 28 i 29 i 30 Here we explore the consequences of 
such asymmetry in detail with the goal of extracting in- 
formation about impurities in graphene from transport. 

The starting point for the calculations is the tight- 
binding model for clean graphen o 17 i 18 
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where t is the nearest neighbor hopping amplitude. The 
operators aj and bj creat electrons on sites i and j in the 
graphene honeycomb lattice. The fact that this lattice 
can be considered as two displaced triangular lattices, 
with a unit cell consisting of two atoms, here denoted A 
and B, is made explicit by the introduction of the two 
creation operators. In reciprocal space, the Fermi surface 
is reduced to two inequivalent so-called K-points at the 
corners of the first Brillouin zone, which we find at K„ = 
Airukx/ (3v3a), where a is the nearest neighbor distance, 
v = ±1, and k x is unit vector along the /c^-axis. At 
low energies e <g; t, the dispersion is linear around these 
points, e(k) ~ ±u/|/c|, where k is the k- vector measured 
relative to the K-point, and v/ — iat/2 is the Fermi 
velocity. The clean limit retarded Green's function in 
the K-point (low-energy) approximation is the inverse of 
the 2x2 (in the space of the atoms A and B) Dirac 
Hamiltonian matrix and has the form 2 ^ 
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where [3 = a,vg{n x +in y ) is an angle defining the direction 
of the vector k with respect to the fc^-axis. 
We include impurities by adding 
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to the Hamiltonian, where V lmp is the impurity strength. 
The number of impurities N lmp in the two sub-lattices A 
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FIG. 1: (a) The impurity self energy and (b) the density of 
states for large impurity potential V\ mp . An impurity band of 
width 7 ~ e c ^/n\ mp is formed near the Fermi level at an energy 
e r ~ tc/Mmp. For £>7 the spectrum is linear. 



and B is assumed approximately equal and small com- 
pared with the number N of unit cells in the crys- 
tal. In the dilute limit, when we only keep terms of 
first order in the density n\ mp = iVj mp /iV, crossing dia- 
grams are neglected when the configuration average is 
performed over the random distribution of impurities. 32 
The resulting self-energy is written in terms of a t-matrix 
± R (e) = n- lmp f R (e), where 



f R (e) = V m 



(V ]mp /N)"£6 R (k,e) 
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The sum over k can be performed analytically in 
the K-point approximation. The off-diagonal compo- 
nents vanish, the resulting self-energy is diagonal, and 
the average Green's function is (1/N) ^ k <j (k, e) = 
l(z/e 2 c ) log [-z 2 /(e 2 c - z 2 )] , where z = e - S R (e). The 
energy cut-off e c is related to the cut-off k c in reciprocal 
space below which the dispersion is linear. 

We plot the self-consistent self-energy and the result- 
ing density of states in Fig. [1] In the unitary limit, de- 
fined as V\ mp — > oo, a band of impurity states is formed 
centered at the Dirac point ^^^i^iao w idth, 
estimated by 7 = — 3E fl (0), is computed by solving 



the equation (27 2 /e 2 ) ln(e c /7) 



To logarithmic 



accuracy, 7 scales as y/n; m p. For V\ mp deviating from 
the unitary limit, the impurity band is shifted away 
from the Fermi level to an energy e r . The self-energy 
is approximated as S (e) = — ij + a(e — e r ) + ... for 
small energies e — e r -C 7. To lowest order, we have 
e r = -e 2 /(2Vi m p[ln(e c /7) - 1]) and a = a /(l + a ), 
where ao = (l/ni mp )(27 2 /e 2 )[ln(e c /7) — 1] is only weakly 
(logarithmically) dependent on the density of impurities. 
The shift of the impurity band away from the Fermi level 
by the amount e r leads to a large electron-hole asymme- 
try and the anomalous thermoelectric response that we 
study below. For large energies, e > 7, the density of 
states is essentially the same as in the clean limit. 

Nearest neighbor hopping it') give an electron-hole 
asymmetric contribution to the dispersion, but it en- 
ters in second order in the low-energy expansion: e(/t) = 



±v f \n\ + v 2 f \n\ 2 [t' /t 2 T f cos(3/3)/(6t)] + ... These contri- 
butions are only important far from the Fermi level and 
we neglect them. 

The linear response is defined as 32 



J 

JE 



£-11 £12 

C21 £-22 



E 



(4) 



where E is the electric field and T is the temperature. 
The Onsager relation £12 = £21 holds. The charge 
conductivity is defined as a = e 2 £n/T, while the elec- 
tronic contribution to the open-circuit heat conductivity 
is K el = (1/T)(£ 22 — £\ 2 /£n). Finally, the thermopower 
is defined as S = — £12/ (eTCu). The response functions 
£ij are defined in terms of Kubo formulas, i.e. current- 
current correlation functions. The charge current opera- 
tor for graphene modeled as above has the simple form 

where h.c. denotes hermitian 
ic. Note that the function d>u — 



j = J2k Vka^k + h.c 
conjugate, and Vk = Vk 



and therefore also Vk, are complex func- 
tions since the three nearest neighbor vectors Sj point 
in three directions rotated 120° relative to each other 
in the graphene plane. The heat current operator has 
an extra term from impurity scattering, j^; = + j'™ p , 
where j° = (l/2)£ k (v k k + KM( a W + b l b *)> and 
j£ P (l/2)E 



kq 



which 



(C/^Vk + £^Vk+q)a k 6 k+q + h.c. 

in principle complicates the calculation of the thermal 
response. However, for graphene, as was shown for nor- 
mal metals by Jonson and Maha n 33 i 34 (see also Refs.l35[). 
the response function kernels are simply related to each 
other: once we know the charge current response ker- 
nel K\\ — K(c) we get the other two kernels through 
K 12 = cK(e) and K22 = e2 K(e). The repsonse functions 
are then computed by integration 
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where /(e) is the Fermi distribution function. For point 
scatterers the kernel K (e) can be approximated by the 
bare bubble ^ 2 ' 23 The needed sum over k can be com- 
puted analytically in the K-point approximation. We 
write z = e - Z R (e) = a(e) + 16(c) and find 2 ^ 
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In Fig. [2] we present results for conductivities and the 
thermopower. The conductivities are only weakly de- 
pendent on the exact value of V\ mp as long as V\ mp e c , 
as is clear when we compare the solid and dashed lines 
in Fig. [2ja)-(b). These response functions are given by 
the electron-hole symmetric part of K(e) which is ro- 
bust as long as the impurity band is not shifted far from 
the Dirac point. We conclude that, for these response 
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functions, the unitary limit is effectively approached 
quickly for Vj mp 3> e c . On the other hand, the ther- 
mopower is given by the electron-hole asymmetric part 
of K(e) and is very sensitive to impurity scattering. We 
have S(— V\ mp ) = — S(V, mp ). The order of magnitude 
S ~ fcs/e ~ 100/iV/K is in agreement with experiments 
on single- wall carbon nanotubes^ 

At low temperatures both the charge conductiv- 
ity and the slope of the thermal conductivity reaches 
constant universal values a = 4e 2 /(irh) and n = 
AirkgT/ (3h), respectively, that are independent of the 
details of the impurities, in agreement with results in 
the literaturej 19 ' 20 i 21 ' 22 i 23 On the other hand, the low- 
temperature slope of the thermopower reaches a non- 
universal constant value, that depends sensitively on the 
nature of impurity scattering. These results are under- 
stood in terms of a Sommerfeld expansions of the trans- 
port coefhcitents in the small parameter T/7 <C 1 for 
slowly varying kernels K{j (e) = (Q) +e (e) | £ _ Q + . . . . 
The needed parts of the kernel Eq. © can be computed 
analytically, K(0) = 2, and K'(0) = 4a(l - a)e r /(37 2 ), 
where we assumed that t r « 7 « e c . The universal con- 
ductivities immediately follows, and the low-temperature 
thermopower is (reinstating h and ks) 

^ 2a(l - a)e r k%T 

37 1 ~- {) 

This means that at low temperatures, ksT <C 7, the 
Wiedemann- Franz law is obeyed, L = n el /aT = Lq — 
(7T 2 /3)(fc 2 3 /e 2 ), and also the Mott formula holds, S = 
~{n 2 /3)(k 2 B T/e) d[lnK(e)}/de\^ Q . At higher tempera- 
tures T > 7, however, neither of these relations hold, see 
Fig. |2Jd)-(e). This happens for any system where the 
conductivity kernel is varying around the Fermi level on 
some particular energy scale, here given by the impurity 
band width 7. For normal metals this scale is typically 
given by the much larger Fermi energy. 

There are other contributions to the heat conductance 
besides the electronic. In particular phonons are impor- 
tant in graphite and carbon nanotubesj 3 "^ for which the 
Lorenz ration is typically found to be L ~ lOLo — lOOLo. 

In contrast, the thermopower is given entirely by elec- 
tronic contributions. In the low-temperature limit, T <C 
7, we see that S is proportional to e r and 7~ 2 , where 
for large V\ mp > e c , e r oc e 2 /Vi mp and 7 cx e cy /n~. 
The thermopower is therefore propotional to the in- 
verses 1/VJmp and l/n.j m p. In the strict unitary limit, 
Vj mp — > co, when electron-hole symmetry is restored, the 
thermopower vanishes i 21 ' 31 But for a large (but not in- 
finitely large) impurity potential the thermopower is en- 
hanced: the smallness of l/VJ mp is compensated by the 
large l/n mp - This also means that if n\ mp can be con- 
trolled, Vi mp can be extracted by measuring the slope of 
the thermopower at low temperatures. 

We note that there are clear analogies with the sit- 
uation in a d-wave superconductor^ where the Lorenz 
ratio deviates from L but recovers at low temperatures, 
T <C 7. An anomalously large thermoelectric coefficient 
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FIG. 2: Temperature dependence of (a) the charge conduc- 
tance, (b) the heat conductance, and (c) the thermopower for 
three values of the impurity density n lmp . Solid lines are for a 
large impurity potential Vi mp = 20e c and the dashed lines for 
the strict unitarly limit Vimp — » 00. In (a)-(b) we have normal- 
ized the conductivities by the low-temperature asymptotics 
0-0 = 4e 2 /(-irh) and k (T) = 4irk%T/(3h). In (d) we show the 
details at low temperatures. The dash-dotted lines are the 
thermopower computed throught the Mott relation. In (e) 
we show the Lorenz ration L = k/((tT) in units of the value 
Lq — (tv 2 /3)(k% /e 2 ) appearing in Wiedeman- Franz law. We 
used e c = leV — * 11605-/^ to convert the temperature scale 
to Kelvin. 



£12 has been predicted? 3 "^ A significant difference for 
graphene compared with the superconducting case is the 
possibility to measure the thermoelectric response di- 
rectly through the thermopower, which is not possible 
in a superconductor since supercurrents short-circuit the 
thermoelectric voltage. 

The chemical potential of graphene can be tuned by 
applying a voltage to the substrate^ We show the de- 
pendence on the chemical potential in Fig. GHa)-(b). The 
charge conductance as function of /i is essentially linear 
at large /i. The value of the thermpower at /1 = and 
the associated asymmetry of S around \x = 0, and also 
around the point //*, where S(fi*) = 0, is related to the 
amount of electron-hole asymmetry caused by impurity 
scattering. 
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FIG. 3: The charge conductance (a) and the thermopower (b) 
as function of chemical potential at T = 50-rY In (c) we show 
how thermal smearing leads to an approximate Mott's rela- 
tion based on the charge conductance as function of chemical 
potential. For T C 7, the Mott relation is exact to leading 
order (dash-dotted line). 



Mott's formula, derived mathematically through the 



Sommerfeld expansion, is only valid at low temper- 
atures T <C 7, since the conductivity kernel K(e) 
varies slowly only on the energy scale 7. But as we 
show in Fig. [21c), it turns out that thermal smearing 
leads to an effective, approximate Mott relation S ~ 
-(tt 2 /3)(fc|T/e)d[lno-(ju; T)]/d/j, in terms of the full tem- 
perature dependent conductivity. The non-linear tem- 
perature dependence S(T) is obtained at /i = 0. 

In summary, we have presented results for the linear re- 
sponse to electric and thermal forces in graphene for the 
case of strong impurity scattering, near the unitary limit. 
The impurity band of width 7 is centered near the Fermi 
level, at an energy e r oc e 2 c /V\ mp . The induced electron- 
hole asymmetry gives small changes of the charge and 
thermal conductivities, but leads to an enhanced ther- 
mopower, which at low temperatures T <C 7 is linear in 
temperature with a slope proportional to l/(Vf mp 7ii mp ). 
The thermopower, measured directly or estimated by 
Mott's relation, can therefore be used to extract infor- 
mation about impurity scattering in graphene. 
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